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Abstract.  Various  methods  for  treating  nonhomogeneous  Dirichlet  boundary  con¬ 
ditions  for  the  p-version  of  the  finite  element  method  are  presented.  These 
methods  are  theoretically  and  comcutationally  analyzed.  Numerical  experimen¬ 
tations  are  given.  They  clearly  illustrate  the  importance  of  the  right  treat¬ 
ment  of  the  nonhomogeneous  Dirichlet  boundary  conditions. 


1 .  I ntroduct ion. 


The  classical  finite  element  method  approximates  the  solution  of  the 
problem  by  piecewise  polynomials  (usually  of  degree  p  =  1,2).  Accuracy  is 
then  achieved  by  decreasing  the  elements  size  h,  which  is  why  this  method  is 
often  called  the  h-version.  The  h-version  has  been  thoroughly  investigated 
both  theoretically  and  computationally  and  many  codes  based  on  this  approach 
are  available. 

In  recent  years,  attention  has  been  focused  on  two  new  finite  element 
developments,  the  p-verslon  when  the  mesh  is  fixed  with  accuracy  being 
achieved  by  increasing  the  degree  p  of  the  elements  either  uniformly  or 
selectively,  and  the  h-p-verslon  which  Increases  the  degree  of  elements  and 
modifies  the  mesh  simultaneously.  The  first  theoretical  papers  addressing 
these  new  versions  appeared  in  1981  ([1)  and  [2])  and  discussed  basic  approxi¬ 
mation  and  convergence  results  for  these  metho'ds.  (See  also  [3],  [4]  for  the 
state  of  art  of  these  methods. )  These  methods  are  related  to  the  spectral 
methods  which  are  now  used  in  computational  fluid  mechanics  (see  eg.  [5]). 

In  structural  mechanics  the  equations  to  be  solved  are  of  elliptic  type 
and  the  input  data  are  piecewise  analytic.  This  piecewise  analyticity  has 
profound  impact  on  the  regularity  of  the  solution.  The  solution  is  piecewise 
analytic  with  singular  behavior  in  a  priori  known  areas  such  as  neighborhoods 
of  corners,  see  [6].  It  can  be  shown  ((7),  [8])  that  for  this  type  of  proDiem 
the  p-versior  *as  a  rate  of  convergence  (in  the  energy  norm)  which  is  twice 
that  of  the  h-version  (with  respect  to  the  number  of  the  degrees  of  freedom) 
when  a  quasiuniform  mesh  is  used.  The  h-p-version  then  leads  to  »n  exponen¬ 
tial  rate  of  convergence.  For  more,  see  [9),  [10]. 

The  p  and  h-p  versions  have  been  implemented  in  the  commercial  code 
PROBE  (Noetic  Technology,  Inc. ,  St.  Louis)  which  is  used  in  the  industries. 


The  program  STRIPE  (Aeronautical  Research  Institute,  Sweden)  is  also  based  on 


the  p  and  h-p  version. 

In  [11],  112],  [13],  we  discussed  a  class  of  theoretical  questions 
related  to  the  implementation  of  essential  boundary  conditions  in  the  frame¬ 
work  of  the  p  and  h-p  versions.  The  main  aspect  is  to  approximate  the 
nonhomogeneous  boundary  conditions  so  that  they  will  be  in  the  trace  spaces  of 
the  finite  element  subspaces.  This  approximation  can  be  made  in  various  ways. 
We  introduce  in  this  paper  a  one  parametric  family  of  approximations  based  on 
the  expansion  in  Jacobi  polynomials  which  naturally  lead  to  the  employment  of 
the  theory  of  weighted  Sobolev  spaces.  We  analyze  this  family  theoretically 
and  perform  various  comparative  numerical  studies. 

Section  2  of  the  paper  introduces  the  weighted  spaces,  and  interpolation 
results  between  families  of  these  spaces.  It  also  contains  a  family  of  pro¬ 
jection  operators  based  on  Jacobi  polynomials  which  are  analyzed  both  theore¬ 
tically  and  computationally  in  section  3.  These  results  are  applied  to  the 
p-version  of  the  finite  element  for  solving  problems  with  nonhomogeneous 
Dirichlet  boundary  data  in  two  dimensions.  The  results  of  numerical  experi¬ 
ments  are  also  presented  in  this  section. 
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2.1.  Notation. 


2 

Let  £1  c  R  be  a  polygonal  domain  with  straight  or  curved  sides.  We 
M  _ 

assume  that  3£1  =  T  =  Ur  where  f  ,  i  =  1,2 .  M  is  an  analytic  simple 

i  =  l  1  1 

arc  connecting  the  vertices  A^_^  and  A^  (Aq  =  A^).  By  we  denote 

f.  -  (A.  .  v  A. )  and  by  w,  the  interior  angle  of  A,,  0  <  w  <  2 n.  The 
l  l-l  i  i  i  i 

notation  scheme  is  shown  in  Figure  2.1.  Further  we  denote  I  =  (-1,1)  c  R*. 


Figure  2.1.  Scheme  of  the  domain  and  the  notation. 


2.2.  The  spaces. 

We  will  consider  the  usual  Sobolev  spaces  LgCG)  =  H^(n),  H^( f2) ,  Hk(fi), 
k  >  0  integer.  The  norm  and  seminorm  will  be  denoted  by  IMI  .  and 

H*(n) 

respectively.  We  will  also  consider  the  spaces  with  k  fractional. 

The  definition  then  is  the  usual  one  based  on  the  K-method,  see  eg.  [14], 
[15]. 


3 


On  I  we  will  consider  a  family  of  weighted  spaces.  For  s  £  0  an 


integer  and  p, v  real  we  define  the  norm 


W^p.tO 


(2.  la) 


+  1 

!u||2 

w(m,u)  J_1 


[(1-x2)  ^f^rl  +  (1-x2)  yu2]dx. 


For  0  <  s  =  [s]  + {s}  where  [s]  is  an  integer  and  0  <  {s>  <  1  we  define 


(2. lb) 


W S(p,i>) 


1 

,^L 


2. -p/2  (s),  ,  ,,  2, -p/2  (s),  .,2 

[(1-x  )  ^  u  (x)-(l-y  )  ^  u  ( y ) ] 


-1J-1 


Ix-yl 


l+2ts) 


(s)  d  u  _  n  j, 

where  u  =  - .  For  s  =  0  we  use  u  =  v  and 

dxS 


|ul20  -  \  I 

VT(p.u)  J_x 


(l-x2)'Mu2dx 


(instead  2j(l-x2)  *\i2dx) 


Let  further  C°°,  respectively  C*.  denote  the  set  of  functions  with  all 
derivatives  on  I  =  [-1,1],  respectively  functions  with  all  derivatives  and 
compact  support  in  I. 

We  define  now  Ws(p,y)  and  Vs as  the  completion  of  the  set 

(u  €  C°°  |  ]|  uj|  <  »> ,  resp.  of  the  set  C*.  The  spaces  W^p.r)  are 

WS(p.u)  U 

Hilbert  spaces.  Ws(0,0)  is  the  usual  Sobolev  space  with  fractional  deriva¬ 
tives,  i.e.,  Ws(0, 0)  =  H®( I ) .  We  will  consider  in  the  sequel  various  spaces 

and  their  equivalency.  The  equivalency  will  be  denoted  by  =.  The  weighted 
spaces  Ws(p,i;)  have  been  studied  in  a  general  context  in  [14].  Let  us 

mention  now  seme  of  their  properties  which  will  be  used  later. 

Theorem  2.1. 


(a)  Let  <s}  *  1/2.  Then  all  spaces  Vs (u,  i>)  with  -®  <  v  <  2s  + 
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then 


are  equivalent .  If  in  addition  p  +  2s  *  1  +  2k,  k  =  0 . [s  ]  -  1 , 

these  spaces  are  equal  to  Ws(^i,p+2s). 

(b)  If  v  >  fi  +  2s  then  W^ji,  *i+2s)  =  Ws(^,p+2s). 

(c)  If  s  =  0  then  ft^(n,i>)  =  W^(ji,v)  for  any  -®  <  v  <  00.  . 

For  the  proof,  see  Theorem  3.2.6  and  Remark  6-3. 2. 6  of  [14].  □ 


Let  (H  H_)  be  the  interpolated  space  by  the  K-method.  For  more 

1  c  6,  c 

see  [14],  [15].  Then  it  is  well  known  that  for  0  <  {s}  <  1,  s  i  0, 


(2.2a) 


HIsMs)  «  <HIs1.HIsM) 


(s) , 2’ 


and  for  {s>  *  1/2 
(2.2b) 


g[s]+{s}  _  (g[s]  o[s]+l} 


{s> , 2‘ 

By  H3  we  denoted  the  standard  Sobolev  spaces  HS(fl)  or  HS(I). 
then 


We  have 


respectively 


HIsMs>(I)  «  WIsl*(s,[0,0). 


Htsl*<s>(l)  =  StS,*<S>(0.0) 


In  the  case  {s>  =  1/2  the  interpolated  norm  H' S  * ' I i  has  the  form 

=  [,',iijj[s]*l/2  *  11  ( >-x2)'1/2“<S)  i  j]  • 


1/2 

This  norm  Is  usually  denoted  by  (I).  We  will  consider  only  the  interpo¬ 
lated  norms  in  this  paper.  Hence  +  0)  not  equivalent  with 

HIs1*I/2(I). 

In  the  sequel  we  will  write  (•,•)  instead  of  (•,•) 

o  6,  ^ 

It  is  possible  to  interpolate  between  the  spaces  W^t^i  ,y  ),  i  =  1,2, 
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and  VIs1  Ip.  .  v  ) ,  i  =  1,2. 


Theorem  2.2.  Let  s ^  2.  0,  s2  £  0,  *  s^t  „  £  ji  +25^  i  =  1,2,  0  <  6  <  1, 


a)  (Mr^)s2  =  ^2~v2)sv 

b)  s  =  ( l-0)s^ + 9s2,  s  *  an  integer, 

c)  v  =  ( 1-0 )  i>  i  +  Qv 2> 

v^i  ^2~v2 

d)  — —  -  — - —  =  -  (in  the  case  that  s,  =  0  and  s„  >  C,  one 

S  S1  S2  12 

H-v  M2~v2 

sets  n=v.  and  uses  -  =  — — -). 

11  s  s„ 


(2.3) 


^2^2'V2])e  = 


For  the  proof,  see  Theorem  3.4.2  of  (14). 


We  will  now  apply  Theorems  2.1  and  2.2  in  special  cases  we  will  be  intej — 
ested  in  later.  Consider  i  =  1,2,  with 


Then  by  Theorem  2. 1 


s,  =  0,  s„=r,  r  >  1  an  integer 


^1  =  Vl  =  V’  0  <  v  <  1 


M2  =  v  -  r,  v2  =  v. 


W°(Mr»i)  *  W°(41,i^1) 


\T {lL?,v0)  =  \ir ( Mg ’ y2 ^  ^ (n0,v0) 


where 


Mi  =  Mi’  ui  =  Mi  +  2si’  1  =  1>2- 
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Further  we  have 


0  =  (VU1)S2  =  (M2_1;2)S1 

and  for  0  <  0  <  1.  0r  *  an  integer 

s  =  ( 1-0^  +  0s2  =  0r 

y  =  ( 1  -0 ) u ^  +  6v2  =  2s  +  [  (  1-0)m1  +  =  s  +  u 

p  =  y  -  2s  =  y  -  s. 

For  is}  *  1/2  and  p  +  2s  =  y  ♦  s  =  v  *  1  +  2k,  k  =  0 . [s]-l  we  have 

fl^p.y)  =  VIs  ( i.  s) 

*  ( W°(fij ,  Uj  j ,  Wr(M2,y2))Q 

•  urC»z,v2))e 
=  (W°(y,y),  W^y-r.y)) 

o 

and  so 

V^Cp.y)  =  (W(y,y),  W^y-r.yUg 

for  y  <  y  =  0r  +  y . 

Assume  now  that  0  <  v  <  1,  v  *  1/2,  r  -  1,  and  0  =  y.  Then  we  have 
s  =  0,  p  =  0  and  hence 

(&°(y,y),  ftX(i;-l,w))w  =  (0, 2y )  =  tf'iO.O)  =  f(I). 

We  see  that  by  interpolating  the  weighted  Sobolev  spaces  W^Cp^.y^),  i  = 
1,2,  we  can  obtain  the  standard  fractional  Sobolev  spaces. 
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2.3.  Polynomial  bases  on  I. 

Let  us  create  (eg.  by  Gramm  Schmidt  procedure)  the  orthogonal  polynomial 
basis  in  the  space  V°(a,a).  Denote  this  set  of  polynomials  by  {P^xia)} 
where  n  is  the  degree  of  the  polynomial. 

Theorem  2.3.  Let  a  <  1.  Then  P^txia)  are  the  Jacobi  polynomials  with  the 
index  /3  =  -a,  i.e.,  P^x;  a)  =  P^x;  -0,  -0)  where  P^x-.-p.-p)  is  the 
standard  Jacobi  polynomial. 

The  theorem  follows  immediately  from  the  basis  properties  of  Jacobi  poly 
nomials,  see  eg.  [lb],  [17],  D 


We  will  list  the  major  properties  of  the  Jacobi  polynomials 


[2.4a) 


Pn(-i.«)  -  <-i>nrn(ii«)  -  l-Dn  ^TT^TT 


(2.4b) 


,1 


(1-x2)  'V  (x; a)P  (x;a)dx  =  0  for  n  *  m, 
n  m 


-1 


(2.4c) 


, ,  2,-«,„  ,  ,,2,  2"2a+1r2(-a+n+l)  _  ,2,^ 

(1-x)  (Pn(x;a))  dx  =  n ;  ( -2«+2n+  !  )  1~  ( -2a+n+T7  ~  Va) 


(2. 4d) 


P'(x;a)  =  ^(-2a+n+l)P  (x;a-l) 

n  2  n- 1 


(2. 4e)  (l-x2)P"(x;a)  -  2(-a+l)xP' (x;a)  +  n(-2a+n+l)P  (x;a)  =  0 

n  n  n 

or 

((l-x2ra+1P'  (x;a))'  +  n(-2a+n-*-l)(l-x2)"aP  (x;a)  =  0. 
n  n 

We  will  also  define  the  orthonormal  Jacobi  polynomials  P^(x;a) 

(2.5a)  P  (x;o)  =  P  (x;a)A-1(a) 

n  n  n 


and  get 


(2.5b) 


P  ( + 1 ; a)  = 
n 


-1)^  (-1,  a)  = 
n 


d  (a)n 
n 


-a* (1/2) 


where  d  (a),  n  =  0,1,2,...  is  bounded  from  above  and  below  by  constants 
depending  only  on  a. 

So  far  we  assumed  that  a  <  1.  Let  us  address  now  the  case  a  =  1 . 
Obviously  a  polynomial  P(x)  belongs  to  W{ 1 . 1 )  only  when  P(±l)  =  0. 

Theorem  2.4.  Let  a  =  1.  Then  the  orthogonal  polynomial  basis  < P  ( x ; 1 ) }  in 

W(  1,1)  is  given  by 


(2.6) 


Pn(x;l) 


r* 

P  , ( t ; 0 )dt ,  n  £  2 
n-l 

-1 


and  (2.4b),  (2.4c)  hold  for  n  >  2. 

Proof .  P^( t ; 0 )  is  the  Legendre  polynomial.  Because  of  the  orthogonality  of 

Legendre  polynomials  ?n f ± 1 ; 13  *  0,  n  >  2.  Obviously  Pn(x;l)  belongs  to 
W°  ( 1 , 1 ) .  We  have  to  show  that 


,1 


(1-x2)  lP  (x; 1)P  (x; l)dx  =  0 
n  m 


-1 


for  n, m  >  1 ,  n  *  m.  For  a  =  0  we  have  from  (2,4e) 


and  hence 


P  (x;  0)  =  -(  ( l-x2)P'  (x;  0) ) '  .  ■ 

n  n  (n+l)n 


P(x;l)  =  -i-(l-x2)P'  (x; 0) . 

n  2n  n-l 


Therefore 

-1 


(l-x2)_1P  (x;l)P  (x; 1 )dx  =  C(m,n) 
n  m 


(l-x2)P'  ( x; 0 ) P'  . ( x; 0 ) dx 
n-l  m-1 


-1 
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=  C^m,  n) 


1 


(1-X  )P  _ ( x ; - 1 ) P  (x; -1 )dx 
n-2  m-2 


when  using  (2.4b,d)  for  m  *  n,  n,m  2  2.  It  is  easy  to  check  that  (2.4c) 
holds  for  a  =  1  and  n  2  2. 


Using  now  (2.4c)  and  (2.4d)  we  see  that  the  system  (Pn(x;a)>  is  an 

1  o  l 

orthogonal  base  in  W(a-l,a)  for  a<l  and  of  W  (a-l,a)  for  a  =  1. 


2.4.  The  approximation  properties  of  the  polynomials  on  I. 

Let  P(x; a) ,  a  £  1,  be  the  set  of  the  orthonormal  polynomials  in 
W^(a,a).  Assuming  u  €  \P(a.,a),  a  £  1,  we  can  write 


(2.7) 


-  ZcA(x;o) 


and  the  series  (2.7)  converges  in  W  (a,a). 


Let  us  define 


(2.8) 


c,  P,  ( x ;  a) 
k  k 


(2.9) 


*u  =  u  -  f au  =  )  c,  P.  (x;a). 
3  p  L  ^  k 


We  have 


i  Pk  ( x;  a )  *1  W°  ( <x- 1 ,  <x- 1 )  =  (a)||Pk(x;a)|lW°(a-l,a-l) 
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so  that 


=  A^*  (a)  A^ia-1  )^-( -2a+k+l )  =  B^(a)*k 


0  <  C^(a)  <  B^(a)  S  02(a)  <  ® 


P^(x;a)  =  (x; a-1 )kD^(a) 


where  D^(a)  is  bounded  from  above  and  below  by  the  constant  dependent  only 
on  a. 

Assume  now  that  u  €  W*(a-l,a)  then  u'  e  W°(a-l,a-l).  For  u  given  by 
(2.7)  we  have  now 


(2. 10) 


u'  *  XckkVk-l(x;“‘,) 


and  hence 


(2.  11a) 


00 

:  (a)  Hull2  <  Yc2(k2+1)  <  c  (a)||u||2 

VT  (a-1 ,  a)  K  Vr  (a-1 ,  a) 


Similarly  for  any  integer  r  >  0  and  u  €  W1" (a-r,a)  given  by  (2.7)  we  have 


(2.11b)  C  (a, r)  ||u|!  <  >  c, 

^(a-r.a) 


i  Y cj(k2r.l )  <  C2(a,  rljul!2 


\F  (a-r,  a) 


(2. lib)  leads  immediately  to 

Theorem  2.5.  Let  r,s  be  integers  and  r  >  s  £  0,  u  €  \f  (a-r, a) ,  a  <  1. 


II  <fu|| 


^  W^a-s.a) 


<  C(a,r,s)p  S ^ || u || 


v/” (a-r,  a) . 


So  far  we  have  assumed  u  €  W1'  (a-r,  a),  s  £  r  and  we  have 


11 


fau  e  W^a-s.a).  Let  us  now  consider  the  case  when  u  €  (a-r,a).  Then 

P 


(X  o  1 

obviously  in  general  ¥  u  t  W  (a-l,a). 


Lemma  2.6.  Let  n  >  0  an  integer,  a  <  1,  and 


I  n/2] 


(2. 12a) 


•  Z  <2k*l)P2k(x;«) 


k=0 

[(n-l)/2] 


(2. 12b) 


4“i(x)  = 

2,  n 


y'.  (2k+l)P2k+1(x;a). 


k=0 


[a] 


Then  ip  (x)  have  following  properties: 
l ,  n 


[a] 


a)  ip  (x)  is  an  even  function  and  ip-  (x)  is  an  odd  function; 


1,  n 


2,  n 


[a], 


b)  \pz  (x),  i  =  1,2  is  a  polynomial  of  degree  n; 
i ,  n 


,  ,  _  .  -a+(5/2)  ^  ^  -a-s(5/2)  ,  ,  „ 

c)  (2.13a)  c.  (a)n  i  ip,  (1)  £  c_(a)n  ,  i  =  1,2 

l  l ,  n  2 


d)  (2.13b)  d.  (a)n(3/2)"'r  <  S  d_(a)n 

i  i , n  ,  2 

w  (a-r,a) 


(3/2)+r 


r  =  0,1,  i  =  1,2. 


The  above  statements  follow  immediately  from  the  basic  properties  of  Jacobi 
polynomials. 


Let  us  now  define 


, [a] ,  . 
ip.  (x) 


C(*%)(lM-l)i  +  1(y"u)(-l))4i£ 

P  P  2  P  P  ..laJ 


i=l 


* 

*  •  p 


and 


2  [oc]  (x) 

*  p'u-iV (<?>%m)*(-i)i,’1(j’au)(-i))-L4 — 

i=l  P  P  *[“}(!) 


Then  obviously  (£au)(±l)  = 

P 


^°au)(±l)  =  0  and  £au,  £eau  €  W  (a-r.a),  r 
P  P  P 
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0,1.  Further,  let  §au  =  u  -  £au  and  §°au  =  u  -  £  au. 

P  P  P  P 

Theorem  2.7.  Let  0  <  a  <  1 ,  r  =  0,  1 ,  s  an  integer,  sir,  a  >  1 


(2. 14) 


||5“u||  £  C(a,r,s)p 

P  W^(a-r.a) 


-(s-r) 


^(a-s,  a) 


Proof.  Because  of  the  definition  of  {^(a-s.a)  it  is  sufficient  to 
(2.14)  for  u  €  c"  Let  s  i.  1,  1  >  a  >  1-s 

00 

"  *  ZcA(x;ai 


and 


k=0 


u( ±1 )  =  0. 


Hence 


(2. 15) 


I 

k=p+l 


|(f*u)(±l)i  S  2.  lckMPk(±l;a)| 


<  C  ^  |cklk_a+(1/2)  <  C  £  |cklkSk 
k=p+ 1  k=p+ 1 


s,  -a+(  1/2)-: 


<  C 


I 

Hc=p+1 


.  ,2.2s 

ICk1  k 


1/2 


-Or+l-S 


-  ^Ilull  _  P 

w(a-s,  a) 


-a+l-s 


Therefore,  using  Theorem  2.5  and  Lemma  2.6  we  get 


il$“u||  <  ||Q“u|| 

p  w(a-r.a)  P  V^Ca-r.a) 


.  r'M  H  -a+l-s  1  (3/2)+r 

*C|lullus,  ,p  -«^5/2>p 

VT  ( a-s, a)  p 


s .  Then 


prove 
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<  Cp-(s~r)  ||u|| 


^(a-s.a). 


Hence,  with  r  *  1, 
>1 


|  [(l-x2)'(a"1)[(6“u)']2-»-  (l-x2)_a(Q“u)2Jdx  5  Cp_(s_1) 


Hull 


VIs  (a-s,  a) 


and  ($pU)(±l)  *  0.  Using  Theorem  330  of  [19]  we  get  for  a  >  0 


,♦1 


-1 


(l-x2r(1'*a)($au)2dx  £  C  (l-x2)“(a”l,l(6au)']2dx 
P  |  P 

' -1 


which  shows  that  also 


ll§“u||  S  C||5“u||  £  Cp”(s_1 5  ||u|| 

P  Wia-l.ct)  P  W  (a-l,a+l)  w(a-s,a) 


and  (2.14)  is  proven  for  r  *  1,  ski. 


Next  let  u  €  W°(a, a)  n  c",  a  <  1.  Then 


ICPpU)  (±1)| 


r 

‘I 


|ckl  |Pk(±l;ot)  | 


r 

“I 


k=0 


|cklk"a+(1/2)  <  CHull  0  p  a+1 

w  (a, a) 


and  hence 


i  Cm 


■„-»  0  ^  IIQ“ull  +  C||u||  Ip"8*1  ■  — TsTpT  P3/2I  *  CHull 

P  W°(a,a)  P  W°(a, a)  W°(a,a)  Pa+(5/2)  W°(a,a) 


Because  yau  €  ft^O.a)  for  a 
P 

1  _  j>1 

P  “  P‘ 


*  1  the  theorem  also  holds  for  a  *  1  when 


Theorem  2.7,  together  with  Theorems  2.1  and  2.2,  leads  to  a  series  of 
important  approximation  results.  As  we  have  seen  earlier 


^(ji.D)  =  (ft°(a,  a) ,  ^(a-r,  a) ) 
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where 


s  =  6r,  s  *  integer,  {s}  *  1/2 
p  -  a  -  Sr 

|i  +  2s  =  v  =  s  +  a*l  +  2k,  k  =  0 . [s]-l 


v  £  v. 


Hence  by  Theorem  2.7  we  get  for  0  <  0j  <  1,  0j  *  1/2,  s  i  1,  an  integer, 
a  £  1,  a  >  1  -  s,  (r  =  1) 


(2. 16) 


iffu 


p"a9i<^ 


S  Cp'(s"Slllu| 


vV 


{^(a-s, a) 


where 


\  *  °'8i 


\  m  “*ei 


and  (2. 16)  holds  (only)  under  the  assumption  that  0^  *  1/2,  s  >  1  an 
integer. 

Let  us  remark  that  we  can  replace  on  the  left  hand  side  of  (2. 16)  the 

004  004 

norm  of  W  l(a-0^,a+0^)  by  the  norm  of  the  space  W  J(a-0^,a).  Let  now  0  < 

0  <  1.  0  *  1/2 

s0  =  Is0]  +0. 


If  [s]  £  1  then  we  can  use  the  interpolation  on  the  right  hand  side  and  get 


l$fu| 


P  "  ,®,9 


W  Ma-Oi.a) 


-(vV 

5  cP  °  iiuii 


(a-sQ,a) 


where  we  have  excluded  the  case 
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a  *  se  =  1  +  2k*  k  =  0 . 1  se 1  "  1  • 


If  [s]  =  0  then  first  by  simultaneous  interpolation  on  both  sides  of  (2. 14) 
we  prove  that 


n$“un  *  cimi 

P  W  1(a-6i,a)  W  1(a-01,a) 


and  then  using  (2.16)  for  s  *  1  by  interpolating  on  the  right  hand  side  we 
get 


l«u| 


W01(a-0i,a) 


£  Cp"(e"0l)||u| 


W0(a-0, a) 


provided  that  0  2  @1  and  0,01  *  1/2.  We  have 

Theorem  2.8.  Let  0  <  a  5  1,  0  <  £  1,  s2  >  s^,  s^  *  [s. ]  +  1/2,  a^  +  s^  * 

1  +  2k,  k  -  0 . [Sjl  -  1.  Then 


(2.17) 


l£ul 

P  WSl(a-sl,a) 


S  Cp*(S2"Sl)||u|| 


^(a-sa.a) 


In  addition  we  can  replace  the  spaces  W^ta-s^.a)  by  W3*  (a-s^  ,  a+s^  )  in 
(2.17)  and  also  WSl(a-s^,a)  by  ^'(g.y)  where  fS  S  a-s^,  y  <  a  and 
WS2(a-s2>a)  by  WS2(3,y)  with  8  >  a  “  s2  ?  >  a.  □ 


(Because  of  direct  use  of  (2. 14)  we  do  not  exclude  the  case  of  s.  being 
integers). 

Select  now  a  =  Sj,  0  <  a  <  1,  a  *  1/2.  Then  we  get  from  (2. 17) 


(2. 18a) 

Select  now 

(2. 18b) 


i$!u« 


p  aa(o,o)  p  sa( i ) 


S  Cp'(Sz‘a)||u| 


ftS2(a-s2,  a) 


st  <  a,  s1  *  1/2  in  (2. 17): 


,l$DU|,0S  -  =  ll^Uli 

P  HSl(I)  P  WSl(0,0)  P  ^'(a-Si.a) 
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£  Cp~(s2"Sl)||u|| 


WS2(a-s2,a) 


For  s^  >  a  we  obviously  cannot  use  the  argument  above.  In  (2.18)  we 

excluded  the  cases  s^  *  [s]  ♦  1/2  and  a  +  s^  *  1  +  2k,  k=  1 . Is^J  —  1. 

In  the  next  section  we  will  be  especially  interested  in  approximation  in 
« 1/2 

the  space  H  (I)  which  was  excluded  from  our  consideration.  Hence  we  will 
consider  and  get  from  (2.18a)  the  following. 

Theorem  2.9.  Let  5  i  e  >  0,  a  =  i+c,  s  >  a.  Then 

tL  c. 


(2. 19) 


>  Y1/2)+cci) 


.  _  -(s-(l/2)-c)..  .. 
*  Cp  u 


Vs  (a-s,  a) 


for  a  +  s  *  1  +  2k,  k  =  0 . (s)  -  1.  □ 

On  the  left  hand  side  of  (2.19),  we  can  replace  fi(1/2)+C(I)  by  H^2(I) 
and  on  the  right  hand  side  ^(a-s.a)  by  f^O.y),  0  2.  a-s,  y  i  a. 

Remark  2.1.  Using  Theorem  2.9  we  have  lost  pe  in  the  estimate  of  the  error 
1/2 

i-n  Hqq  (I).  This  case  can  be  studied  separately.  For  example  in  [11]  we 
have  shown  that  for  s  >  1/2 


||$°1/2u|| 

D  1/2 

P 


.  r  -(s-(  1/2) ) ,  .1/2.,  , 

^  Cp  (log  p)  l|u|| 


^(0,0) 


So  far  we  have  addressed  the  approximation  properties  of  Let  us  now 

analyze  §  a.  We  have 
P 

2  r-Ia),  )  .[“If  *•. 

rau-^“u  =  \ y  ((?“uhi)  +  (-i)l+1(f“uHD)  4^ — 14—  • 

p  p  h  p  p 


By  (2.  15),  we  have  for  s  >  -a  +  1  21s  before 


I  (7%)  (±1)|  $  Cp"a+1"S||u| 


Vs  ( a-s ,  a ) 
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Remembering  that  (£°au  -  ^u)  (±1 )  =  0  using  Theorem  2.1  and  also  Theorem 
330  of  [19]  we  get  for  a  >  0 


°au-£au)||  *  C||(Kau-£“u)|| 

p  p  VT(a-l.a)  P  P  W1  (a-1 , a+1 ) 


C||rau-^u)|| 

p  p  &1  ( a- 1 ,  a ) 


Hence  by  interpolation  we  get  for  0  S  i  1, 


_  <;  cp-a+1-S2[upa'(5/2)p(3/2)+Sl]||u||oi; 

P  P  WSl(a-s1,a)  WS2(a-s2,a) 


k°au  -  *au|| 


£  Cp“(S2_Sl)|tu|| 


£fS2(a-s2,a) 


provided  -a+1  5  ana  as  before  a  +  *  1  +  2k,  k  *  0,  .  .  .  ,  [s2]  -  1. 

Let  vis  summarize  our  results  in 


Theorem  2.  10.  Let  0  <  a  £  1,  Sj  1  0,  Sj  +  a  *  1  +  2k,  k  =  0, . 
i  =  1,2.  Then 


.As)-  1, 


1)  for  s2  >  Sj,  Sj  *  lsi]  +  ^ 


(2.20a) 


l$“u|  S  Cp_(s2_Sl)  Hull 

P  WSl(a-Si,a)  WSz(a-s2,a) 


If  in  addition  -a+1  <  s^  then 


(2.20b) 


||3°au||  <  Cp'(s2_Sl)  Hull 

P  WSl(a-Si,a)  WSz(a-s2,a) 


2)  For  0  S  0  S  a  S  1,  s2  >  0,  0  *  1/2 


(2.20c) 


l$Ju|  S  Cp"(S2“P)||u|| 

P  H^(I)  WSz(a-s2,  a) 


If  in  addition  -a+1  £  0  <  s2*  then 
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(2.  20d) 


||^au||  <  Cp"(  S2‘0 }  ||u|| 

P  H*(I)  fi^U-sz.a) 


It  was  essential  in  (2.20c)  (2.20d)  that  3  £  u.  The  case  a  <  0  has  to  be 
studied  separately.  Let  us  mention  first 

Lemma  2.11.  Let  u  be  a  polynomial  of  degree  p  and  0  £  a  <  3-  Then 

||u||  S  Cp2(P_a)||u|!  .  c 

ir(i)  h  (i) 

For  the  proof,  see  eg.  Lemma  2.4  of  [18]  or  ll). 

Now  we  can  prove 

Theorem  2.12.  Let  0  <  «  S  p  S  1,  #,M  1/2,  3  <  s2>  s2  *  ts  ]  +  s2  +  a  * 

1  +  2k,  k  =  0 . ts2l  -  1.  Then 


(2.21a) 


ll$“u||  <  Cp'(s2"0)  +  (P'a)||u|| 

P  fl*(I)  f2((3-s2,{i) 


If  in  addition  1  -  3  <  a  <  s2  then 

(2.21b)  ||§°au|l  <  c  -(s2-3)  +  (3-a) 

P  H0(I) 


^O-Sz.p) 


Proof.  Using  (2.20c)  we  get 


s  cp-(S2-a) 


p  fid) 


lull 


VfS2(a-s2,  a) 


l4?ui  s  cp-(s=-a,i. 

P  fid)  fiS2(g-s2,  S) 


and  hence 


Itti-^UI  =||^au-^u||  <  Cp"(S2"a)||u|| 

P  P  fta(I)  P  P  Ha(I)  ^(p-sz.p) 


Using  Lemma  2.11  we  get 
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and  hence 


llS“u-$?u||  ft  5  Cp2(0~a)  (s2"a)Hu|| 
p  p  H^(I)  WS20-s2,  P) 


,,3“ui  ,cp-(S2-0)+(^a)iiu|iOiS9 

P  fl^d)  V^20-s2.P) 


where  we 


wrote  on  the  left  hand  side  ( I )  instead  of  H^(I)  because  3  * 


1/2,  u  €  H^(I)  and  u  is  a  polynomial  vanishing  at  x  =  ±1.  The 

P  P  P 

second  part  of  the  theorem  follows  easily,  too.  c 


Remark  2.2.  Using  Theorem  2.12,  we  see  that  for  a  >  1/2,  s  £  (1/2)  +2e, 
c  >  0,  we  have 


«<SUV/2 

Hoo 


C(c)p-<s-‘1/2))*clul 


6s (a-s, a) 


o  o(X 

and  the  same  estimate  holds  also  for  Q  • 

P 

In  the  next  section  we  will  especially  be  interested  in  the  approximation 
1/2 

in  the  space  HqQ  ,  which  case  was  excluded  from  our  consideration.  We 
addressed  this  case  only  via  the  approximation  in  C(I). 

It  is  well  known  that  IP^(x; 1/2)}  are  Tchebyshev  polynomials  which  can 
be  written  in  the  form 


P^(x;l/2)  =  cos(k  arc  cos  x) 


Pk(x;l/2) 


•  /I 


cos  (k  arc  cos  x) . 


We  have  shown  in  [11]  that  if 


u  *  Zw— 


k=0 


then 
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00 


(2.22) 


1/2 

H17  1 1 ) 


=  1 


lckl2k*  lc0l2. 


k=0 


1/2 

The  space  H^g  (I)  then  has  a  norm  which  is  equivalent  to 


(2.23) 


hoo2(i) 


■1 


(1-X2)”1u2dx+  Hull2 

H1  ^(1) 


This  norm  is  equivalent  to 


I  ,  .  Obviously  *1/2u  and  £°1/2u  € 

61/2(o, i)  p  p 


Hg22(I)  If  u  €  ftk(I).  k  >  1/2. 


In  [11],  we  have  proven  various  properties  of  £°*/2  for  example 


n  A°  1/2  I.  „  -  -(k-(  1/2) ) ,  1/2  „ 

II  Q_  u||  ,  <  Cp  (log  p)  ||u|| 

P  Hqq2  ^ I )  Hk(I) 


k  >  1/2. 


Let  us  remark  that  in  general 


,/2 

P  *oo"] 


and  ll<3p1/2u|| 

P  H^2(l) 


(and  analogously  in  other  cases)  are  not  necessarily  monotonic  non- increasing 

functions  of  p.  Of  course,  if  we  define  as  the  projection  onto  the 

set  of  polynomials  of  degree  <  p  in  the  norm  (2.23),  then  the  monotonicity 

1/2 

would  be  guaranteed.  If  we  neglect  the  first  term  in  (2.23)  and  define  f 


as  a  projection  in  this  norm,  then  we  get  a  monotone  sequence  but  now 

,1/ 

P 


?1/2  *  H1/2 

*  Hoo  • 


2.5.  Numerical  experimentation. 

In  the  next  section  we  will  analyze  the  error  of  the  finite  element 

method  for  two  dimensional  elliptic  problems.  This  analysis  leads  naturally 

to  the  measure  of  the  error  on  I  in  the  norm  II  •  II  _  ,  respectively 

H1  (I) 
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Hence  we  will  be  interested  here  in  the  computational  analysis  of 


"  "  1/2 

H00  (I) 

the  error  in  the  one-dimensional  case,  in  the  norm 


1/2 

H  (I) 


More  precise¬ 


ly,  we  will  analyze  II u  -f  u  || 

P  H  (I) 


where 


u  (x)  = 
w 


'l-2wx+x  k=0 


P^(x;  0)  ,  0  <  Cl)  <  1. 


We  choose  the  functions  u^(x)  since  they  are  characteristic  for  applica¬ 
tions.  In  practice,  in  the  finite  element  method,  the  singularity  is  almost 

always  located  at  an  end  point  of  I.  If  u — »1  then  u  (x)  becomes  singu- 

u> 

lar  at  x  =  -1.  Hence  the  parameter  w  characterizes  the  smoothness  of  u 

u> 

when  the  nonsmooth  behavior  occurs  at  an  end  point  of  I. 


Ilu-A.!  ... 
r>>  -  -  ■’  "  »> 

'V'J  1/2 

P  H  (I) 


where  II  •  II  ,  is  given  by  (2.22).  Obviously  ra(w)  >  1  and  ra(u) 

TT  *  '  ^  ✓  T  ^  D  D 


H  (I) 

expresses  the  quality  of  the  performance  of  ?a  in 

P 


1/2 

n  7  l  I  ) 


Figure  2.2 


shows  r^fw)  for  u  =  0.9,  0.95  and  0.98  and  for  the  (extreme)  cases 

a  =  0,1.  Table  2.1  gives  the  values  r°(0.9)  and  ^(0.9).  The  values  of 

P  P 

Tp(0.95)  and  r^(0.98)  are  close  together. 


Table  2.  I. 

The  values  ra(0.9),  a  =  0,  1 
P 


p 

r°(0.9) 

P 

r^O.S) 

P 

1 

1 . 0000 

1 . 0000 

2 

1.0113 

1.1611 

3 

1. 0904 

1.2741 

4 

1. 1067 

1.2450 

5 

1. 1516 

1.3218 

10 

1.2602 

1.3281 

15 

1.3192 

1.3854 

20 

1.3892 

1 . 3720 

25 

1. 4078 

1. 4038 

30 

1. 4663 

1 . 3909 

35 

1.4676 

1.4130 

40 

1 . 5207 

1.4018 

50 

1. 5622 

1.4089 

We  see  that  r^(a)  decreases  (in  cur  range)  as  a — >1  while  r\a) 

P  P 

increases.  For  p  large  (with  respect  to  r-~)  the  projection  P^  performs 

1-(j>  p 

better  than  P^,  but  for  practical  p  (p  <  15,  say),  P^  seems  to  be 

preferable  to  P*.  As  we  will  see  in  the  next  section  a  similar  effect  is 
P 


present  in  the  context  of  the  finite  element  method. 
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3.  The  Finite  Element  Method 


3.1.  The  model  problem. 

Let  £2  be  the  domain  as  given  in  Section  2.1.  Let  M  =  {1.....M}  and 

V  c  M.  Denote  now  r0  =  U  f  and  =  T  -  r°.  We  will  call  rD,  respect- 

i€® 

ively  r*  the  Dirichlet,  respectively  the  Neumann  boundary.  Let  H  =  H  (£1) 
and  HQ  =  {u  e  Hlu  =  0  on  rD>  and  B(u.v)  be  a  continuous  symmetric 
bilinear  form  on  HxH.  We  will  assume  that  for  all  u  e  HQ  we  have 

B(u,u)  £  rl|u||jj,  7  >  0. 

Further,  let  F  be  a  continuous  linear  functional  on  H  and  g  e 

g[l1  =  g|  €  Hk(T  ),  i  €  2),  k  >  1/2. 
i 

Our  model  problem  is  now 


(3.  1) 


’Find  uQ  €  H  such  that 
_D 

uQ  =  g  on  T 

and 

B ( Uq , v )  =  F(v) ,  V  v  €  Hq 


Our  model  problem  has  a  unique  solution. 


3.2.  The  p-version  of  the  finite  element  method. 

q  _ 

Let  £2  =  U  £2  ,  where  £2  ,  i  =  1 . q  are  open  curved  quadrilaterals 

i=l  1 

or  triangles  called  elements  of  the  partition  of  £1.  The  vertices  of  £L  are 

called  the  nodes  of  the  partition.  We  will  assume  that  the  nodes  which  are 

located  on  the  boundary  T  of  £1  coincide  with  the  vertices  A.  of  £1  as 

J 

shown  in  Figure  3.1. 
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A4  y 


Figure  3. 1.  The  scheme  of  the  partitioned  domain. 


By  S  =  (-1,1)  and  T  -  {£,tj|-1  <  £  £  0,  0<t)<  ($+1)v5;  0  <  £  <  1, 

0  <  i}  <  ( 1  — C ) v^}  (see  Figure  3.2)  we  denote  the  standard  square  and  standard 

triangle  respectively.  Further,  let  P*(T)  be  the  set  of  all  polynomials  of 

2 

degree  £  p  on  T  and  P^IS)  the  minimal  set  of  polynomials  on  S  consist 

ing  of  all  polynomials  of  degree  5  p  and  polynomials  which  are  of  degree 

5  p  on  one  side  and  are  zero  on  the  three  others.  For  details  see  eg.  (20). 
2 

The  set  Pp(S)  is  the  set  of  serendipity  elements  and  in  [21]  is  denoted  by 


84  I  1  B| 

I  *  I  ■  I 


(o,o)  \  e 


Figure  3.2.  The  scheme  of  S  and  T 


Assume  now  that  the  mappings  :=  (xj^  =  xj*^  (€,*!).  i  =  1,2),  J  = 

l,...,q  are  smooth  one-to-one  mappings  of  S  on  £2.  if  £2.  is  a  (curved) 

J  J 

quadri lateral  and  T  on  £2j  if  £2^  is  a  (curved)  triangle.  We  can  now 
speak  about  vertices  and  sides  of  £2^  in  the  obvious  sense. 

We  will  assume  the  following  about  the  partition  of  £2  and  the  mappings 


i)  The  intersection  £2^  n  £2^  is  either  empty  or  is  the  single  common 

vertex  of  £1  and  n,  or  is  the  single  entire  side  of  £2  and  £1,. 

ii)  If  -  fij  n  Qj  is  the  common  side  of  £2j  and  £2^,  P  c  P  = 

F,  (pj)  *  Fj(Pj),  Pi  €  BkBk+1 •  pj  e  ,  (Bj  is  the  vertex  of  S  or  T), 

then  d(P^,B^)  =  dCP^.B^)  where  we  denoted  by  d(P^,B^)  the  Euclidean  dis¬ 
tance  between  P.  and  B,  . 

i  k 

We  also  will  identify  the  side  with  I  *  (-1,1)  and  the  map  = 

*  • 

Fj  =  Fjj  of  I  onto  riy  Realizing  that  the  sides  of  T  and  S  have 

•  • 

length  =  2  the  relation  between  F^  F^,  F^  is  obvious. 

Let  now  P  (£1)  =  (u  €  H1  (£2)  |  u( F. 1  (x, ,  x_ ) )  6  P1  ( T)  if  £1.  is  a  (curved) 

p  1  1  2  p  1 

-1  2 

triangle,  u(F^  (x^.x^))  €  Pp(S)  Oj  is  a  (curved)  quadrilateral).  If 

u  e  Pp(£2)  then  u(Fi 1 (x^, x^) )  is  a  polynomial  of  degree  p  on  every  side  of 
T,  respectively  S.  We  will  identify  the  sides  of  T,  respectively  S, 


with  I  =  (-1,1)  in  the  obvious  way. 

Assume  now  that  a  mapping  is  given  which  maps  H1^  I),  k  >  1/2,  onto 

Pp(I),  (i.e.,  the  set  of  polynomials  of  degree  p  on  I),  with  ($pU)(±l)  = 

u(±l)  and  &  u  =  u  for  u  e  P  (I).  We  note  that  for  k  >  1/2, 

P  P 


H^d) 


c — »C(I)  and  hence  u(±l)  is  well  defined. 


The  p-version  of  the  finite  element  method  for  solving  problem  (3.1) 


consists  of  finding  u  €  P  (£2)  such  that 

P  P 
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(3.2a) 


(3.2b) 

where 


B(u  , v)  =  F(v)  V  v  e  P  (Q)  o  Hn, 
P  P  0 


up  *  $pi]glil  on  rr  i  €  X, 


*!lVi](x)  *  F*$  (g[iI(F*_1(x))) 
P  i  P  1 


and  F^  is  the  mapping  of  I  on  induced  by  the  mapping  F^. 

The  following  theorem  gives  an  error  estimate  for  the  p-version. 

Theorem  3,1.  Let  the  solution  u^  of  the  problem  (3.1)  have  the  form 


(3.3a) 


(3.3b) 


M  nJ 

V1  V  til  ti] 

uo  =  ui*L  Lcj  uj 

j=i  i=i 

u1  €  I ^(Q).  k  >  1, 


(3.3c)  uj»  *  rJi^Uog  rJ|3rJ113Cjil(rJ)*jilCOJ) 


where  (r^.G^)  are  the  polar  coordinates  with  the  origin  at  A^,  a^i+1^  > 


atlJ  >0, 

s) 


.  Z^(r,)  is  a  C°°  cut-off  function  and  is  a  C°° 


function  in  0  .  Then 


(3.4) 


|U0'UdI!  1  *  c[p'(k_1)||u  ||  +p-fi|log  p|S  Y  |c[i 

P  IT(Q)  t-  1  H*(n)  f-  J 

1  »  J 

ts>  p  *4(ri)J 


where  p  *  min(2oJ.^),  s  =  where  p  =  2a\^\  C  is  a  constant  indepen- 

i.j  J  0  J 

dent  of  Uq  and  p  but  dependent  on  k,  p,  s  and  on  the  partition  of  Q. 
Proof.  The  proof  of  (3.4)  is  essentially  the  same  as  of  Theorem  4.1  in  [11]. 
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Let  us  mention  that 


H1/2{r  ) 

Hoo  (V 


(I) 


where  g^*(£)  =  £  €  I,  because  by  the  assumption  F^C)  is 

smooth. 

Theorem  3.1  shows  that  the  error  of  the  finite  element  solution  depends 
on  two  terms.  The  first  two  terms  in  (3.4)  represent  the  error  of  the  best 
approximation  without  the  effect  of  imposing  the  Dirichlet  boundary  condition. 
In  the  case  when  D  *  <p  or  g^  =  0  only  these  terms  are  present  in  (3.4). 
The  last  term  shows  the  effect  of  the  approximation  of  the  Dirichlet  condition 
in  the  dependence  on  the  projection  &  .  We  will  call  this  term  the  con¬ 
straint  term.  We  defined  in  Section  2.4  the  projections  and  £°a  for 

the  set  of  functions  u  which  satisfy  the  condition  u(±l)  =  0.  In  order  to 

choose  &  *  £a  we  define 

P  P 

a<x  •  -*  Aa,  *  -V 

*  u  =  u  +  r  (u  -u  ) 

P  P 

*  *  _* 

where  u  (£)  =  u(F^(£)),  £  e  I,  and  u  is  the  linear  interpolant  (in  £) 

*  *  _*  AO<x 
of  u  so  that  (u  -u  )(±1)  =  0.  The  mapping  is  defined  similarly. 

For  a  simplification  of  the  notation  we  will  write  often  £a  and  $,°a 

P  P 

instead  of  $a  and  $°a  . 

P  P 

We  can  now  directly  combine  Theorems  2.9,  2.10,  2.12,  and  Theorem  3.1  to 

get  the  error  estimate  for  the  p~verslon  of  the  finite  element  method  when 

using  $a  *  or  &°a  «  9>°a  .  For  example,  let  1  £  a  >  1/2,  s  >  o,  $°a  * 

P  P  P  P  P 

o.a 

.  Then  the  constraint  term  has  the  estimate 
P 
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where  ^  and  g^ ^  were  defined  above.  We  assume  that  s  *  Is]  1/2 

and  a  +  s  *  1  ♦  2k,  k  =  0 . [s]  -  1.  We  obviously  “lost"  a  power  pG  in 

(3.5)  in  comparison  with  the  optimal  estimate. 

We  have  assumed  that  a  >  1/2.  If  a  <  1/2  then  we  have  to  apply 


Theorem  2.12  and  get  on  the  right  hand  side  of  (3.5)  the  term 

—  ( s—  ( l/2)-e)+(  ( 1/2) -a) ..  [i]*  -til* 
p  llg  _g  IL 

W( (l/2)+e-s, (l/2)+c) 


3.3.  Numerical  experiments. 

Let  us  consider  the  case  of  an  L-shaped  domain  shown  In  Figure  3. 1.  We 

will  partition  this  domain  Into  three  squares  as  Indicated  also  in  Figure  3.3. 

du  du 

a4  r4/dn  a3  r*/3"  a2 


U  /  a, 

u 


Ae  /  r7  A 1 


■du  FOR  NEUMANN 
dn  PROBLEM 


Figure  3.3.  The  L-shaped  domain. 


Let  us  consider  the  model  problem 


(3.6a) 


Au  =  0, 


(3.6b) 


u  =  g  on  T1  and  fg  (i.e.,  1,8  €  2?), 


(3.6c) 


aH  =  h  on  rr  2*1*7. 
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We  will  assume  that  g  and  h  are  such  that  the  exact  solution  u^  is 
(3.7)  Uq  =  r*  sin  A0 


where  ( r,9 )  are  the  polar  coordinates  with  the  origin  at  the  vertex 
(see  Figure  3.1).  Note  that  u^  *  0  on  This  problem  will  be  called  the 

Dlrlchlet  problem.  In  addition  we  will  consider  the  Neumann  problem  when 
D  -  {1>.  The  problem  now  is  formulated  in  the  form  (3.1)  with 


B(  u,  v) 


* 

F(v)  =  hvdx  i  €  M.  -  2). 

Fi 

The  constraint  term  is  not  present  in  the  Neumann  problem.  It  is  present  on 
Tg  in  the  case  of  the  Dirichlet  problem 

We  will  now  use  the  projection  which  is  advantageous  in  practice 

(although  in  some  cases  some  assumptions  in  our  theorems  Eire  violated). 

The  Neumann  problem  has  no  constraint  term  and  hence  the  error  for  the 
Dirichlet  problem  has  to  be  larger.  The  difference  indicates  the  influence 
of  the  constraint  term  (and  of  on  the  accuracy  of  the  solution.  In 

Tables  3.  1  -  3.5  we  show  the  relative  error  |u_-u  I  ,  /|u_|  .  for 

0  p  h’isi)  0  h'iq) 

various  a  and  A. 

We  see  that  in  the  case  A  =  1.6  and  A  =  0.  6  the  ratio  of  the  error  of 

the  Neumann  problem  and  the  Dirichlet  problem  is  a  reasonable,  one  for  all  a 

and  nearly  independent  of  p.  In  [7]  we  have  proven  that  in  this  case,  when 

[8]  1 

8  €  H  (Tg),  the  rate  for  the  Neumann  and  Dirichlet  problem  is  the  optimal 

one  when  ^  (or  was  used.  We  cannot  theoretically  explain  why  the 
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error  decreases  with  a — >0  for  a  <  1/2.  Very  likely  it  relates  to  the  used 

range  p  <  8  and  analogous  effects  we  have  seen  in  Figure  2.2. 

[  8  ]  1 

In  the  case  X  -  0.  15  when  g  t  H  (Tg)  we  see  the  divergence  when 

a  o  1  Aol/4 

"p  is  used,  as  expected,  but  once  more  we  cannot  explain  why  y  gives 

o ©0  a  ©  1  /A 

best  results.  We  see  that  t  performs  slightly  worse  than  f  (but 

P  P 

still  better  than  £°*^). 

P 

©  o0 

In  the  case  X  =  0.05,  the  use  of  r  leads  to  the  divergence. 

A  natural  question  arises.  What  a  should  be  used  in  practice?  The 
experiment  suggests  the  choices  a  5  1/2.  a  =  1/2  should  be  preferred  for 
theoretical  reasons  but  a  *  0  usually  gives  the  best  results  in  our  experi¬ 
ments  (these  and  others)  and  never  gives  a  bad  result.  We  note  that  for  X  = 
0.05,  the  p-version  gives  unacceptable  results  even  for  the  Neumann  problem 
and  hence  the  h-p  version  has  to  be  used. 

If  the  h-p  version  (or  p-version  with  strongly  refined  meshes)  is 
used,  then  the  difference  between  various  projection  operators  is  not  too 

important.  We  will  illustrate  one  such  choice  of  operator  when  Dirichlet  con- 

1/3 

ditions  are  imposed  on  the  entire  boundary.  The  solution  is  u  =  r  sin  0/3 


2/3 

(case  A)  and  u  =  r  cos  20/3  (case  B).  The  projection  is  used. 

Figure  3.4  (case  A)  and  3.5  (case  B)  show  the  relative  error  for  different  p 

and  the  meshes  which  are  strongly  refined.  They  have  n  layers  and  for  n  = 

1,  the  mesh  shown  in  Figure  3.2.  Theoretically  (see  [9],  [10]),  the  h-p 

-y?N 

version  converges  exponentially  as  e  *  and  hence  Figures  3.4  and  3.5  are 
plotted  in  the  log  x  scale.  We  see  also  in  case  A  that  the  h-p 

version  with  n  *  p  converges  in  fact  as  e”7  . 


t°l 
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RELATIVE  ERROR  IN  THE  ENERGY  NORM  Hell 


Table  3.1.  Relative  error  for  A  =  1.6 


p 

§>°° 

P 

£°l/4 

P 

£°l/2 

P 

£°3/4 

P 

£°1 

P 

NBC 

1 

2.6001  E-l 

2.6001  E-2 

2.6001  E-l 

2.6001  E-l 

2.6001  E-l 

2.5761  E-l 

2 

3.6361  E-2 

3.6297  E-2 

3.6232  E-2 

3.6171  E-2 

3.6127  E-2 

3. 1326  E-2 

3 

3.0482  E-2 

3.0422  E-2 

3.0367  E-2  ' 

3.0331  E-2 

3.0340  E-2 

2.3647  E-2 

4 

1.5320  E-2 

1. 5361  E-2 

1.5426  E-2 

1.5529  E-2 

1.5699  E-2 

1.2535  E-2 

5 

8.2123  E-3 

8.2512  E-3 

8.3240  E-3 

8.4285  E-3 

8.5922  E-3 

7. 1002  E-3 

6 

4.8521  E-3 

4. 8886  E-3 

4.9435  E-3 

5.0260  E-3 

5.1533  E-3 

4.2937  E-3 

7 

3.0914  E-3 

3. 1174  E-3 

3. 1571  E-3 

| 

3.2182  E-3 

3.3137  E-3 

2.7834  E-3 

8 

2.0821  E-3 

2. 1029  E-3 

2. 1338  E-3 

2. 1809  E-3 

2.2546  E-3 

1.8983  E-3 

Table  3.2.  Relative  error  for  A  =  0.6 


P 

§>°° 

P 

£°l/4 

P 

£°l/2 

P 

£°3/4 

P 

£°1 

P 

NBC 

1 

3.0707  E-l 

3.0707  E-l 

3.0707  E-l 

3.0707  E-l 

3.0707  E-l 

2.6397  E-l 

2 

2. 1588  E-l 

2. 1671  E-l 

2. 1795  E-l  1 

2.2008  E-l 

2.2436  E-l 

1.7145  E-l 

3 

2. 1377  E-l 

2. 1513  E-l 

2. 1720  E-l 

2. 2267  E-l 

2.2781  E-l 

1.5878  E-l 

4 

1.5851  E-l 

1.5991  E-l 

1. 6203  E-l 

1. 6563  E-l 

1.7289  E-l 

1.2927  E-l 

5 

1.2400  E-l 

1.2528  E-l 

1. 2726  E-l 

1.3065  E-l 

1.3763  E-l 

1.0516  E-l 

6 

1.0116  E-l 

1.0233  E-l 

1. 0414  E-l 

1.0725  E-l 

1. 1378  E-l 

8.7797  E-2 

7 

8.5071  E-2 

8.6129  E-2 

8.7780  E-2 

9.0682  E-2 

9.6911  E-2 

7.4923  E-2 

8 

7.3162  E-2 

7.4134  E-2 

7.5652  E-2 

7.8361  E-2 

8.4305  E-2 

6.5076  E-2 
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Table  3.3.  Relative  error  for  A  =  0.3333 


p 

£°° 

P 

£•1/4 

P 

£•1/2 

P 

£•3/4 

P 

£•1 

P 

1 

7. 1888  E-l 

7. 1888  E-l 

7. 1888  E-l 

7. 1888  E-l 

7. 1888  E-l 

5.7124  E-l 

2 

5.8300  E-l 

5.8538  E-l 

5.9181  E-l 

6. 1187  E-l 

7.0125  E-l 

4.7273  E-l 

3 

5.8462  E-l 

5. 9263  E-l 

6. 1026  E-l 

6.5729  E-l 

8.3462  E-l 

4.5779  E-l 

4 

4.9954  E-l 

5. 1145  E-l 

5. 3729  E-l 

6. 0323  E-l 

8.3153  E-l 

4. 1310  E-l 

5 

4.4324  E-l 

4.5536  E-l 

4.8347  E-l 

5. 5783  E-l 

8. 1688  E-l 

3.6981  E-l 

6 

3.9661  E-l 

4.0957  E-l 

4.4009  E-l 

5.2086  E-l 

8.0127  E-l 

3.3549  E-l 

7 

3.6438  E-l 

3.7673  E-l 

4.0747  E-l 

4.9119  E-l 

7.8886  E-l 

3.3079  E-l 

8 

3.3501  E-l 

3.4756  E-l 

3.7924  E-l 

4. 6603  E-l 

7.8025  E-l 

2.8517  E-l 

Table  3.4.  Relative  error  for  A  =  0. 15 


P 

£°° 

P 

£"1/4 

P 

£•1/2 

P 

£•3/4 

P 

£•1 

P 

NBC 

1 

1.0123  E-0 

1.0123  E-0 

1.0123  E-0 

1.0123  E-0 

1.0123  E-0 

8.5288  E-l 

2 

9. 1644  E-l 

9. 1240  E-l 

9.0754  E-l 

9.0775  E-l 

1.0672  E-0 

7.8644  E-l 

3 

9.0381  E-l 

9.0073  E-l 

9. 0034  E-l 

9.2422  E-l 

1.3136  E-0 

7.7035  E-l 

4 

8.4599  E-l 

8. 4380  E-l 

8. 4864  E-l 

8. 9877  E-l 

1.4699  E-0 

7. 3960  E-l 

5 

8. 1341  E-l 

8. 1065  E-l 

8. 1748  E-l 

8.8487  E-l 

1.5893  E-0 

7.0663  E-l 

6 

7.7966  E-l 

7.7722  E-l 

7.8714  E-l 

8.7191  E-l 

1.6836  E-0 

6.7779  E-l 

7 

7.5917  E-l 

7.5626  E-l 

7.6744  E-l 

8.6455  E-l 

1.7626  E-0 

6.5296  E-l 

8 

7.3550  E-l 

7.3276  E-l 

7.4608  E-l 

8.5618  E-l 

1.8308  E-0 

6.3139  E-l 

35 


Table  3.5.  Relative  error  for  X  =  0.05 


p 

£°° 

P 

NBC 

1 

1.0306  E-0 

9.6485  E-l 

2 

i  1.0851  E-0 

9.4059  E-l 

3 

1.2706  E-0 

9.3285  E-l 

4 

1.4422  E-0 

9.2138  E-l 

5 

1.5993  E-0 

9.0854  E-l 

6 

1.7407  E-0 

8.9661  E-l 

7 

1.8707  E-0 

8.8591  E-l 

8 

1.9920  E-0 

8.7630  E-l 
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